arXiv:l504.07312vl [cond-mat.str-el] 28 Apr 2015 


Emergent loop-nodal s±-wave superconductivity in CeCu2Si2: similarities to the iron-based 

superconductors 


Hiroaki Ikeda,'’^’EI Michi-To Suzuki,^ '* and Ryotaro Arita^’^ 

' Department of Physics, Ritsumeikan University, Kusatsu 525-8577, Japan 
^Department of Physics, Kyoto University, Kyoto 606-8502, Japan 
^RIKEN Center for Emergent Matter Science (CEMS), Wako, Saitama 351-0198, Japan 
‘^CCSE, Japan Atomic Energy Agency, 5-1-5 Kashiwanoha, Kashiwa, Chiba 277-8587, Japan 
^Department of Applied Physics, University of Tokyo, Tokyo 113-8656, Japan 
(Dated: April 29, 2015) 

Heavy-fermion superconductors are prime candidates for novel electron-pairing states due to the spin-orbital 
coupled degrees of freedom and electron correlations. Superconductivity in CeCu 2 Si 2 discovered in 1979, 
which is a prototype of unconventional (non-BCS) superconductors in strongly correlated electron systems, still 
remains unsolved. Here we provide the first report of superconductivity based on the advanced first-principles 
theoretical approach. We find that the promising candidate is an s±-wave state with loop-shaped nodes on 
the Fermi surface, different from the widely expected line-nodal d-wave state. The dominant pairing glue is 
magnetic but high-rank octupolar fluctuations. This system shares the importance of multi-orbital degrees of 
freedom with the iron-based superconductors. Our findings reveal not only the long-standing puzzle in this 
material, but also urge us to reconsider the pairing states and mechanisms in all heavy-fermion superconductors. 
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Some sort of intermetallic compounds, containing elements 
with 4f or 5f electrons, are called heavy fermion materials be¬ 
cause the low-energy excitations can be described by heavy 
electrons with the large effective mass up to 1000 times the 
free-electron mass. In these f-electron systems, the strong 
electron correlation and many degrees of freedom bring about 
a variety of curious and interesting physical phenomena, such 
as complex magnetic/multipole ordering, quantum critical 
phenomena, and unconventional superconductivity and so on. 

The heavy-fermion material CeCu 2 Si 2 is the first uncon¬ 
ventional superconductor discovered in 1979 by F. Steglich et 
al. |[T1]. The specific heat coefficient C/T ~ 0.75J/mol-CeK^ 
at Tc ~ 0.5K indicates that the heavy-fermion state has been 
formed by the strong electron correlation between f electrons. 
The BCS-like specific heat jump at is a clear evidence of 
the Cooper pairing of the correlated electrons, in other words, 
the gap opening in heavy-electron bands with high density of 
states (DOS). The relatively high Tc to the effective Fermi en¬ 
ergy Tp ~ lOK made a great impact on the research field 
since it has been considered that the strong Coulomb repul¬ 
sion disturbs the superconductivity. Also, much small Tp, as 
compared with the Debye temperature ~ 200K, is contrary to 
the case in the conventional BCS theory. Since the early stage, 
these facts have implied that the superconductivity is not con¬ 
ventional. Indeed after a short time, unconventional behavior 
in superconducting phase was observed in some experimental 
works; the T-linear behavior at low T in C/T 12 and[2l, and 
no coherence peak just below Tc and the T^ behavior in NMR 
relaxation rate 1 /Ti M. These observations are inconsis¬ 
tent with the exponential decay in the fully-gapped s-wave su¬ 
perconductivity in the BCS theory, rather indicate line nodes 
on the Fermi surface as in the d-wave superconductivity in the 
subsequently discovered high-Tc cuprates. Thus it has been 
widely expected that the line-nodal d-wave state is the promis¬ 
ing candidate of the superconductivity in CeCu 2 Si 2 . Actually 
this material and the high-Tc cuprates share some character¬ 


istic featuresitS; non-Fermi liquid behavior above the optimal 
Tc, and the similar phase diagram in which superconductiv¬ 
ity appears in close proximity to the antiferromagnetic (AFM) 
phase. 

These facts imply also that the Coulomb repulsion, espe¬ 
cially a magnetic fluctuation, is crucially important to make 
mobile electrons bound strongly. In the conventional BCS su¬ 
perconductors, the major glue for electron pairings is quanta 
of lattice vibration, i.e., phonon. The obtained s-wave or¬ 
der parameter has a large weight on the on-site pairing am¬ 
plitude in real space. Therefore, the on-site Coulomb repul¬ 
sion disturbs the on-site pairing and suppresses Tc. On the 
other hand, a plausible pairing glue in the high-Tc supercon¬ 
ductors is the AFM fluctuations driven by the Coulomb repul¬ 
sion, which leads to anisotropic pairs like d-wave state. This 
anisotropic pairing state without on-site amplitude is not di¬ 
rectly suppressed by the on-site Coulomb repulsion. This is 
a reason why the Coulomb repulsion can lead to the high- 
Tc superconductivity. The superconducting pairing function 
is closely related to the pairing mechanism. Clarifying the 
pairing state is crucially important to understand the pairing 
mechanism. Such microscopic studies will give us some use¬ 
ful hints for how to raise the transition temperatures. Indeed 
a recent strategy for finding new high-Tc materials is to in¬ 
crease the energy scale of dynamical spin fluctuations. It is 
based on the fact that the Tc of Cel 15, Pul 15 and the cuprate 
systems is linearly scaled to the characteristic energy of AFM 
fluctuation s 

On the other hand, finding new electron-pairing mechanism 
can offer another route to raise Tc, such as AFM fluctuations 
instead of phonon. In this sense, the finding of the iron-based 
superconductors^^ may give us a chance to search for another 
mechanism, since the orbital degrees of freedom are important 
ingredients in these materials. In this regard, heavy-fermion 
superconductors can provide a playground to search for much 
variety of pairing mechanisms and pairing states. However, 
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the complicated electronic band structure has prevented our 
understanding of nature of superconductivity so far. Under 
the circumstances, recent developments on the hrst-principles 
calculations based on the density-functional theory (DFT) 
have promoted a breakthrough on studying the electronic 
state in the strongly correlated electron systems. For in¬ 
stance, DFTh-DMFT (dynamical mean-field theory) meth¬ 
ods allow us to discuss the correlated electrons, compared 
to angle-resolved photoemmision spectroscopy (ARPES) re¬ 
sults in real materialsi^. On the other hands, as a comple¬ 
mentary method, based on an effective model obtained from 
the DFT calculations, low-energy magnetic/multipole fluctu¬ 
ations have been investigated in the context of the hidden- 
order transition in URu 2 Si 2 Q. The same approach can be a 
powerful tool to dissect the complex electronic state in heavy- 
fermion materials. We here provide the hrst report of the mi¬ 
croscopic analysis about the superconducting gap function in 
CeCu 2 Si 2 . We find that the situation in CeCu 2 Si 2 is similar 
to that in the iron-based superconductors, and the promising 
pairing state is s±-wave state, which is in stark contrast to 
the widely expected d-wave state. The dominant pairing glue 
is magnetic but high-rank octupolar fluctuations. In addition 
we discuss that the second dome of superconductivity under 
high pressuresi^ can possess the similar s-wave pairing, but 
mediated by different mechanism. 

The Fermi surface and magnetic/multipole fluctuations — 
The Fermi surface topology, which is crucially important for 
unconventional superconductivity, can be experimentally clar¬ 
ified by the quantum oscillation and ARPES measurements 
and so on. However, the Eermi surface in CeCu 2 Si 2 is not 
so clear yet. Theoretically, from the electronic band struc¬ 
ture calculations, two possible Eermi surface topology have 
been proposed, which are clearly different each other. One 
is obtained by the ordinary EDA calculations^^, and the other 
by the renormalized band methodic. The Eermi surface in 
the latter is similar to Eig.l(a), which was obtained by our 
LDAh-C/ calculation s'^!*^ . In the EDA Eermi surface shown 
by Eig.l(b), there are a cubic-like electron sheet and a tiny 
electron sheet around F and a complex hole sheet. In the 
LDAh-{ 7 case, there appears a corrugated-cylinder electron 
sheet around X point with heavy effective mass and a com¬ 
plex hole sheet and a tiny hole ring. The band crossing the 
Fermi level is mainly composed of an f-orbital manifold with 
the total angular momentum j = 5/2. In the LDAh-[/ calcula- 
tion^i, it is especially dominated by its azimuthal component 
jz — ±3/2 mixed with a small weight of jz = ±5/2, al¬ 
though the latter components may become larger by including 
the electron correlations22i21. 

First of all, let us investigate the Fermi surface topology 
and magnetic fluctuations. We compute low-energy magnetic 
fluctuations in both EDA and LDAh-C/ Eermi surfaces within 
the random-phase approximation (RPA).^^. In the LDAh-C/ 
Eermi surface shown by Eig.l(c), we find a peak structure at 
around incommensurate Q = (0.21,0.21, 0.5) due to the nest¬ 
ing property in the corrugated heavy-electron sheet around X 
point, which is successfully consistent with a characteristic 
Q vector observed in the neutron scattering measurements^. 
Similar trend has been discussed in a simple model^. On 


FIG. 1. (color online), (a, b) Fermi surface colored by the Fermi 
velocity, (c, d) in-plane magnetic (dipole) RPA susceptibilities for 
<1 ~ {<lx,qy, 0.5), and (e, f) a complete set of multipole susceptibil¬ 
ities along the high-symmetry line^. Left figures correspond to the 
LDAh-C/ case, and right to the ordinary LDA case. In (a), the pres¬ 
ence of heavy-electron sheet around X point is similar to the case 
of the renormalized band methodic. Incommensurate peak positions 
Q — (0.21, 0.21, 0.5) in (c) are consistent with the characteristic Q 
vector observed in neutron scattering measurements^, while in the 
LDA case, (d), it is hard to enhance magnetic fluctuations even for 
larger interactions. In (e) we can see that octupole (Rank 3) fluc¬ 
tuations are remarkably enhanced in the LDAh-C/ case. In the LDA 
case, (f), on the other hand, non-magnetic quadrupole/hexadecapole 
fluctuations are relatively enhanced. Here, susceptibilities without a 
peak at X point are uniformly colored red. 

LDAH-t/ LDA 



the contrary, a hump structure around X in Fig. 1(d) in the 
LDA case is not so much enhanced even for larger interac¬ 
tions. These facts imply that the Fermi surface in CeCu 2 Si 2 
should possess a feature of corrugated-cylinder electron sheet 
as observed in the LDAh-{ 7 and the renormalized band case. 
Hence, we consider hereafter that the Fermi surface at the am¬ 
bient pressure is like the LDAh-C/ Fermi surface hereafter. 

Next let us study a variety of inter-orbital fluctuations^. 
Generally speaking, in multi-orbital systems, there exist the 
orbital degrees of freedom in addition to spin degrees of 
freedom. The presence of strong spin-orbit coupling in 
the f-electron systems entangles these degrees of freedom, 
which leads to the so-called multipole degrees of freedom. 
In the j — hjl multiplet with six jz components {jz = 
±5/2, ±3/2, ±1/2), there exists 6 x 6 = 36 multipole de¬ 
grees of freedom, which are classified into monopole (rank 
0), dipole (rank 1), quadrupole (rank 2), octupole (rank 3), 
hexadecapole (rank 4), and dotriacontapole (rank 5) by sym- 
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metries in the group theory. Conventional charge and mag¬ 
netic degrees of freedom correspond to monopole and mag¬ 
netic dipole, respectively. High-rank multipoles represent 
a spin-orbital coupled degrees of freedom. “Rank” means 
changeable difference of total angular momentum. For exam¬ 
ple, an operator /+ 5 / 2 /- 5 / 2 > which means a change between 
= -1-5/2 and —5/2, is classified into a kind of the highest 
rank 5. This state has been discussed as a promising hidden- 
order parameter in URu 2 Si 2 Q. 

Fig. 1(e) depicts a complete set of multipole susceptibili¬ 
ties in the LDAh-C/ case of CeCu 2 Si 2 . We can see that oc- 
tupole fluctuations are dominantly enhanced, and next is rank 
5 fluctuations. This is consistent with the fact that the ma¬ 
jor component of the ground-state f multiplets is jz — ±3/2 
mixed with small components of= =f 5/2. These high-rank 
fluctuations are much larger than the magnetic dipole fluc¬ 
tuations. This implies that the incommensurate AFM order 
observed in the A-type materials possesses a sizable weight 
of the octupole moment with the same irreducible representa¬ 
tion. On the other hand, in the LDA case shown in Fig. 1(d), it 
is hard to enhance magnetic fluctuations even for larger inter¬ 
actions. Rather, non-magnetic quadrupole/hexadecapole fluc¬ 
tuations become relatively large (Fig. 1(f) )^. Recently, L.V. 
Pourovskii et al. proposed a possibility of orbital transition 
at around the second superconducting dome under high pres¬ 
sures^. In our case this corresponds to a kind of Lifshitz 
transition from the LDAh-C/ Fermi surface with anisotropic 
f-electron charge distribution to the LDA Fermi surface with 
almost isotropic distribution. This implies that with applying 
pressure, incommensurate AFM fluctuations are suppressed, 
but non-magnetic orbital fluctuations are enhanced instead. It 
is an interesting open problem how these non-magnetic fluc¬ 
tuations are related to valence fluctuationsi^ii, which is one 
of hotly-debated issues in heavy-electron systems. 

Superconducting pairing symmetry — Here, let us discuss 
what kind of superconductivity can emerge by multipole fluc¬ 
tuations obtained above. As usual, we evaluate the linearized 
gap equation in the multiorbital systems^, 

^ ^ ^ ^ ^ ) 

k' t't' m'm" (1) 

)Grn'm" {k ) (^k ), 

where /S.im{k) represents the pair function, and Vim,e'm'iq) 
is the pairing interaction. 

Information of multipole fluctuations in this system is con¬ 
tained in xtni,e'm'iq), which is for *e second- 

order perturbation and xf'm^'m' (?) for the RPA. First of all, 
in the RPA susceptibilities, we obtain a dj. 2 _y 2 wave state 
(A = 1.73) and an s± wave (A = 1.43) for the LDA-tf/ Fermi 
surface, which are due to the remarkable octupole fluctuations 
(Fig. 1(e)), while an wave (A = 1.16) and a dxy wave 
(A = 0.56) for the LDA Fermi surface due to the comparably 


enhanced non-magnetic fluctuations (Fig. 1(f)). The leading 
dx‘ 2 -y 2 wave state obtained in the LDAh-( 7 case is the widely 
believed line-nodal d-wave state. However, the sub-leading 
wave also has a large eigenvalue, and it is the leading state 
in the LDA case. In general, within the RPA, magnetic suscep¬ 
tibilities are relatively too much enhanced as compared with 
non-magnetic (charge/orbital) ones. Therefore, the second- 
order perturbation can provide a complementary information. 
Indeed it will be more appropriate in this material, since this 
material has only weak correlation between magnetic fluctua¬ 
tions and transition temperatures. 

In the second-order perturbation, we obtain the s± wave 
symmetry as the leading pairing state, independent of the 
Fermi-surface topology. The leading s± wave has A = 1.06 
and the sub-leading dx 2 -y 2 wave has A = 0.97 for the 
LDA-tf/ case, while the leading pairing state remains the s± 
wave (A = 0.75) for the LDA case. As shown in Fig.2(c), the 
s±-wave state obtained in the LDA-tfA case has a complicated 
structure with loop-shaped nodes on the Fermi surface. The 
similar loop-nodal s± wave state has been discussed in the 
iron-based superconductor BaFe 2 (Asi_a;Pa ;)2 iItHi . Although 
this loop-nodal s±-wave state has still nodal excitations on 
the corrugated heavy-electron sheet, it should be noted that a 
gap size on a flat part is much smaller than that on a convex 
part. With a slight mixture of the on-site pairing due to intrin¬ 
sic attractive on-site interaction by higher-order vertex correc¬ 
tions^, the loop-nodes can be lifted since this nodal feature is 
not symmetry-protected. In this case, the corrugated heavy- 
electron sheet becomes fully-gaped, and only the light-hole 
sheet possesses loop-nodes. This is also the case in the LDA 
case. Thus, the advanced first-principles calculations are in¬ 
dicative that the s±-wave pairing with the loop nodes on the 
light Fermi surface is the promising superconducting state in 
CeCu2Si2. 

Experimentally, recent specific-heat measurements^^ show 
the exponential behavior below 60mK and the i7-linear de¬ 
pendence under the magnetic field H. It has been indicated 
that these features can be explained by multi-gapped s-wave 
state along with the T -linear behavior above 60mK. No co¬ 
herence peak in the NMR 1/Ti does not contradict an s±- 
wave state. Since the peak structure observed in the neutron 
measurements^ is located at around ~ 2 A with the estimated 
gap A, it is not so clear whether this is the resonance peak. 
Rather, such weak peak structure seems to support our high- 
rank multipole-fluctuation mechanism. Thus the s±-wave 
pairing state can be the most probable candidate in this ma¬ 
terial. Then, why such possibility has been missed? First, 
before the discovery of iron-based superconductors, it has 
been considered that power-law behavior in a wide temper¬ 
ature range is a strong evidence of nodal structure, rather than 
multi-gapped s-wave. Second, it has not been so seriously 
taken due to strong sample dependence that its extrapolation 
to OK of a power-law fit in the specific heat is negative. Thus, 
recent improvement of sample quality and low-temperature 
techniques has triggered reconsideration of superconductivity 
in this material. 

Our present results strongly depend on the Fermi-surface 
topology. Therefore, it is crucial to clarify the Fermi surface 
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FIG. 2. (color online), (a) Schematic pressure-temperature phase diagram in CeCu 2 Si 2 d. (b-e) Superconducting gap structures obtained 
in Eq.(l). The bottom figures are the top view or the side view, (b) d^ 2 _y 2 -wdi\e (^cos{2kx) — cos{2ky)) and (c) s±-wave ('^cos{2kx) -I- 
cos(2fcy)) obtained in the LDA-l-(7 Fermi surface, (d) Another s±-wave state in the LDA case, cos{kx) cos{ky) cos{kz), and (e) d^y- 
wave (~sin(fea;) sin(fcy) cos{kz)). Two s±-wave states, (c) and (d), are the promising superconducting state in CeCu 2 Si 2 , different from the 
widely believed d-wave state. The dominant pairing interactions come from the octupole fluctuations in (c) and the quadrupole/hexadecapole 
fluctuations in (d). In addition, these two s±-wave state with different pairing mechanisms may correspond to two distinct superconducting 


phases observed under pressures in (a). 


(a) 

A 



experimentally. In addition, whether the loop nodes exist on 
the hole sheet with light effective mass requires further studies 
at low temperatures in terms of several experiments sensitive 
to low-energy excitations, such as magnetic penetration depth 
and thermal conductivity and so on. 

Finally, we realize a possibility that curious two supercon¬ 
ducting domes (Fig.2(a)) observed under high pressures^ can 
be explained by two types of s±-wave states along with a kind 
of Lifshitz transition from LDAh-[/ Fermi surface to the LDA 
Fermi surface. These pairing states are mediated by two dif¬ 
ferent fluctuations; one is magnetic octupolar fluctuations and 
another is non-magnetic hexadecapolar fluctuations. This im¬ 
plies that the multi-orbital degrees of freedom are crucially 
important. This feature is very similar to two superconducting 
domes recently discussed in the iron-based superconductors^. 
Thus, our study is also helpful to the iron-based superconduc¬ 
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tors. Our findings of the s±-wave state in CeCu 2 Si 2 will stim¬ 
ulate us to reconsider the pairing states and mechanisms in all 
of heavy-fermion superconductors in detail. 
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